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INVESTIGATION OF THE LIFT, CENTER OF PRESSURE, AND 
DRAG OF A PROJECTILE AT A MACH NUMBER OF 8.6 
AND A REYNOLDS NUMBER OF 17 MILLION 


By Thomas N. Canning 


SUMMARY 


Models of a high-speed projectile configuration were tested in free 
flight to determine the lift-curve slope, center of pressure, and drag 
at Mach numbers near 8.6 and at a Reynolds number of 17 million. These 
results were compared with predictions based on available theory and 
experimental results. The drag of a similar cone model was also measured. 


INTRODUCTION 


The aerodynamic consequences of improving the static longitudinal 
stability of a 5° half-angle cone by adding a short spike of heavy mate- 
rial onto the nose are considered in this paper. The resulting shape 
(see fig. 1) was enough different from & cone that some changes in &ero- 
dynamic properties were expected. The lift on the added cylinder was 
expected to increase the lift at angle of attack and to cancel some of 
the gain in stability by moving the center of pressure forward. It was 
also expected that there would be at least a small region of quasi two- 
dimensional flow just behind the spike which might alter the Lift, center 
of pressure, and drag. The presence of the shock wave at the spike base 
was expected to promote boundary~layer transition. The net effects of 
these changes were of sufficient interest that the present experimental 
investigation was undertaken. Tests were made to determine lift, drag, 
and center of pressure of the projectile. The drag of unmodified cones 
was also measured for comparison. 


Since the speeds of interest for this body were in the high super- 
sonic Mach number range, the tests were run at the highest Mach number 
compatible with model strength considerations. The tests were made at 
speeds near M = 8.6 in the Ames supersonic free-flight wind tunnel at a 
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Reynolds number of over 17 million, based on model length. Іп this facil- 
ity the models are fired upstream through the test section of an M = 2.0 
supersonic wind tunnel to attain high relative airspeed. 


SYMBOLS 


a projection in the x,z plane of the instantaneous accelera- 
tion of model center of gravity normal to the tunnel center 
line, ft/sec* 


CD dr&g coefficient, = 
| бао 
Do drag coefficient at а 20 
Or, lift coefficient corresponding to mean vector sum of a and p 
à in test interval 
ас 

Ст lift-curve slope, 4, per radian 
Cm pitching-moment coefficient about model center of gravity, 

pitching moment | І 

Sido 

C dCm 
ma, pitching-moment-curve slope, 2! per radian E: 
D drag, lb 
E,F constants defining variation of œ with time 
e base of Naperian Logarithms 
f frequency of pitching motion, cps 
I . principal moment of inertia about lateral axis through model 

center of gravity, ft-lb sec? | 
k damping constent, вес? 
L | lift, lb 
1 model length, ft 
M test Mach number 


m model mass, slugs 
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Q double integral with respect to time of the angle-of-attack 
function, sec? - 

йо free-stream dynamic pressure, lb/sq ft 

R Reynolds number based on model length and free-stream conditions 

S frontal area of model, sq ft 

t time, sec 

wW projection in x-z plane of the instantaneous velocity of model 


center of gravity normal to tunnel center line, ft/sec 


Хер, Хер distance from model nose to center of gravity апа center of 
pressure, respectively, ft 


х,у; 7 orthogonal coordinate system using tunnel center line ав х 
axis, ft | 

а, angle of attack of model relative to local flight path, radians 

Gg mean value of vector sum of а and В during test interval, 
radians 

В angle of sideslip of model relative to local flight path, 
radians 

e exf, radians / вес 

Subscripts 
1525854 stations in wind tunnel 


APPARATUS AND TECHNIQUES 


Wind Tunnel 


The tests were conducted in the Ames supersonic free-flight wind 
tunnel. This facility is a short ballistic range inside a variable- 
pressure, supersonic, blowdown wind tunnel. While the tunnel was operat- 
ing at a Mach number of 2.0, the models were fired upstream through the 
15-foot-long test section at 5600 feet per second from a gun located in 
the diffuser. The aerodynamic data were obtained from a history of the 
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model motion, as recorded by seven shadowgraph stations (four parallel 

to the x,z plane and three parallel to the x,y plane) and a chrono- 
graph. The data from the four-station group &re used to study the motion 
of the models as if the model motion were confined to the x,z plane. 
Details of the data analysis and a discussion of measures taken to account 
for motion outside the x,z plane are included in the section entitled 
"Data Reduction." Details of tunnel operation are given in reference 1. 


Models 


The models (fig. 1) were machined from 75 ST-6 aluminum bar stock. 
For the projectile the basic 5° half-angle cone was altered by adding a 
cylinder of one third the length, and of 0.19 the base diameter of the 
original cone. The cylinder and nose cone would presumably be made of 
heavy material in a full-size missile to promote stability. The models 
for the present tests were bored out at the base to provide the same 
effect. A conical hole with rounded bottom was bored in each projectile 
model to place the center of gravity at about the Tl-percent point aft of 
the nose. This resulted in a wall thickness of about 0.035 inch, which 
permitted the model to be accelerated about 400,000 gravities in the gun 
without deformation. The nose cones were finished carefully апа the 
bluntest one fired was 0.003 inch flat at the tip; most tips were less 
than 0.001 inch flat. 


The surface finish was produced by a fine finishing cut and a little 
hand sanding in the lathe. No effort was made to improve the surfaces 
further and they were satin-like in appearance. Several model surfaces 
were examined in order to evaluate the equivalent roughness at large 
scale. A photograph of a typical model is given in figure 2, and photo- 
micrographs of a model surface and profile are given in figure 3. The 
pictures are representative of the worst conditions noted; for example, 
the dip, which is indicated in the profile picture is 0.0001 inch deep by 
0.0030 inch long. This profile is believed to be representative of the 
test models, but no extensive examination of model surfaces was made. 

The models were fired from a high performance .50 caliber gun and were 
held on the sabot by evacuating the space in the base of the model as 
indicated in figure 1. The high-pressure powder gases in the gun chamber 
penetrated the plastic-tape seal, and a small amount of gas flowed into 
the evacuated model. The resistance to acceleration held the model firmly 
in place until the sabot was free of the gun. At this moment the powder- 
gas pressure in the model forced separation from the sabot face. | А 
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Data Reduction 


The lift-curve slope, center of pressure, and drag were calculated 
from the record of model motion as a function of time. The component of 
the resultant aerodynamic force parallel to the 2 axis, that is, the 
lift, was obtained by measuring the acceleration along this direction; 
the center of pressure was found from the naturel pitching frequency and 
the liftecurve slope; and the drag was determined by measuring the axial 
deceleration. Details of the lift and center-of-pressure calculations are 
presented below; the details of drag calculation, which is more readily 
visualized, are given in reference 1. For the purpose of data reduction, 
linear aerodynamic characteristics are assumed. The effect of nonline- 
arity is treated approximately. 


As mentioned previously the data from the four shadowgraph stations 
recording the model motion in the х,2 plane were used as if the model 
were free to move in this plane only. The instantaneous acceleration of 
the model projected on the x,z plane may be written 


az _ Chgdo8 


ae a A 


a= 


which is Newton's law of motion in the x,z plane. This assumes a linear 
lift curve and neglects the contribution of the small lift force due to 
pitching and plunging. 


The time variation of angle of attack in equation (1) can be obtained 
using the angle-of-~attack measurements from the shadowgraph pictures and 
two assumptions which define the form of the motion. The assumptions are 
that the restoring moment and damping moments are proportional to the angle 
of attack and pitching rate, respectively. This leads to a variation of 
angle of attack given by 


qe a (B cos wt - F sin wt) (2) 


Equation (2) is fitted to the observed variation of a with respect to 
time by a Least-squares procedure described in reference 2. In this way 
the four unknowns, w, k, E, and F, are evaluated. 


Combining equations (1) and (2) gives the plunging acceleration as 
a function of time 


825 ди _ “ТолоЗ „ку 


一 一 一 = E cos at - F sin ast 
dt? ^ at m (Econ ) (3) 
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Integrating with respect to t gives the following equation for vertical 
velocity, w 


бтеїо8В p* -kt 
те са Wi P aii e (E cos wt - F sin wt) dt (4) 
dt m t 
1 


in which ж is the component of the model velocity along the 2 axis 
&t the first station. А second integration gives the equation of verti- 
cal position relative to the initial point 21. The integrals are evalu- 
ated between limits corresponding to the times and positions in two shad- 
owgraph stations. 


Стадо5 pt2 pte 
2р = 22 + на (te - tı) — | / е "КЕ (E cos ut - F sin wt) аъаъ 
Фа ta (5a) 


The integral on the right can be evaluated since the integrand and limits 
are known. For brevity, it will be designated ©з». 


Z2 = Zi + Wy (te - Бъ) + 12 


Сі,,4о9 | = 

Q (5b) 
In equation (5b) the unknowns are Ст, and иу; the constant wj is elim- 
inated by using data from a third station: 


G S 
23 = 21 + Wi (tg - ti) + 9207 Qi (5c) 


solving equations (5b) and (5c) simultaneously yields the following expres- 
sion for С: 


t2 - t 
m Z-a- (2а - 21) үөт 
Le, = Sq ta - ty (6) 
9 біз - 9з ———— 
tg - tı 


In order to determine the center-of-pressure position, it was also 
necessary to measure the pitching-moment-curve slope, Ста» for use in the 
equation relating lift- and moment~curve slopes: 


= Ung (7) 
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Ihe value Cm, was obtained from the frequency of oscillation ав 
determined by the least-squares fit to the angle-of-attack history. The 
equation for the frequency of oscillation with one degree of freedom was 


used: 
f = 1 = il. -Спа4о51 
= ex? ом l 


w= I 
бш, = = Saot 9) 


The method of calculation of Lift requires data from three stations 
in the wind tunnel to determine one lift answer. For each of the tests 
from which lift-curve slope wes obtained in this investigation, four shad- 
Owgraphs were made. This permitted making four combinations of three 
stations from each set. Occasionally, one combination of three stations 
failed to define the lift-curve slope because the position of the stations 
along the flight path was such as to record three values of z which fell 
nearly on a straight line. Such badly defined values аге not included. 
Typical data from which lift and center-of-pressure location were deter- 
mined are shown in figure 1. 


or 


Thus far the motion has been treated as an oscillation in the x,z 
plane. Actually, the models were free to move in the x,y plane as well. 
lf the aerodynamic properties of the model were perfectly linear, there 
would be no influence on the data if this fact were ignored. Since it 
was expected, in this case, that the lift would be influenced by sideslip 
angle, the value of Cr, was related to the mean absolute angle of yaw, 
the vector sum of а апа В, during the interval in question, rather than 
to а alone. Thus, 


where C is called tbe lift coefficient corresponding to the mean 
absolute angie of yaw, Cig is the value of lift-curve slope obtained ав 


described above, апд ар is the mean absolute angle of yaw for the test 
interval. 


ACCURACY 


A realistic estimate of the accuracy of the above techniques must 
Of necessity be based on consistency of results, rather than on a detailed 
study of all factors contributing inaccuracies. The lift results are 
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consistent within a scatter in Сір of +0.008, based on frontal area, 
for the present test. This amounts to less than i! percent of the meas- 
ured lift at 5° angle of attack. The drag is believed accurate within 
2 percent. The total scatter in center-of-pressure position of six sep- 
arate tests was £1.2 percent of the body length. The static stability 
margin of the models was made small, from 0.3 to 2.5 percent of 1; 
therefore the determination of center-of-pressure position was not sen- 
sitive even to large errors іп (Cm. 


RESULTS 


The testing technique is such that varying amounts of disturbance 
are imparted to the models as they are launched. These differences, plus 
slight damage suffered by a few models, made it necessary to consider 
each test individually. Only the aerodynamic parameters best defined by 
each test have been included here as explained below. 


Lift 


The variation of lift with angle of attack, Орр versus Фр, as deduced 
from four tests (1, 2, 3, and 5), is presented in figure 5. Out of a 
total of 16 lift-curve-slope determinations, 2 were deleted, l each from 
tests 2 and 5, because the data were such as to fail in defining the 
flight-path curvature accurately. 


Test 1 had a peak amplitude of oscillation of about 1° and deviated 
about 0.07 inch from straight-line flight. Tests 2, 3, and 5 had larger 
amplitudes. The model in test 2 was bent about 1.5° at the base of the 
spike during the launch (see fig. 6(b)). The effect of this damage was 
estimated by calculating the theoretical lift on the conical spike nose 
(ref. 3) at 1.5° angle of attack and doubling this to allow for lift _ 
carry-over on the spike cylinder. This lift was then added in coefficient 
form to the values of бір determined from this test and amounted to а 
correction of about 7 percent of the final answer. Test 3 had the largest 
amplitude recorded in the investigation - about 10°. The nose tip (about 
l percent of the model length) was hooked over by lateral acceleration in 
the gun. The hook was bent in a plane nearly normal to the x,z plane, 
so it is believed that this damage had negligible effect on the lift. 


Center of Pressure 


Tests of five models (1, 5, 7, 8, and 9) showed the center of pres- 
sure to be between 72.1 and 74.6 percent of the length aft of the nose. 
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For all of these tests the oscillatory motion was fairly well defined 
and the models were not visibly damaged (see fig. 6). 


Drag 


Five tests were obtained with accurate enough time-distance measure- 
ments for determining drag (1, 4, 5, 6, and 7). These results are plot- 
ted in figure 7 as a function of the average position of boundary-layer 
transition along the model as observed in the shadowgraphs. The drag 
due to lift of these models was estimated and subtracted as outlined. in 
reference 1. 


The principal evidences of turbulent boundary-layer flow are eddies 
visible in shadowgraphs and small waves sent out by these eddies. А region 
of turbulent flow may be seen near the base of the model in figure 6(j). 
This patch of turbulent flow is being swept downstream at about half the 
free-stream velocity as evidenced by the angle of the weak shock wave 
emanating from it. Serutiny of all the pictures from a test permitted 
the estimation of transition position. In the case of test 6 no laminar 
flow was seen aft of the spike cylinder. It is believed that transition 
did, in fact, occur at the nose due to some slight damage not visible in 
the shadowgraphs. Іп view of the difficulty of seeing boundary-layer 
turbulence near the nose, the point for this test is plotted as if tran- 
sition occurred at the spike base. 


The drag tests of two unmodified cones are also reported in this 
figure. The lower Reynolds number of these two tests resulted from using 
shorter models and testing at slightly lower tunnel pressure for testing 
convenience. 


DISCUSSION 


The initial lift-curve slope predicted using reference 3 is indicated 
in figure 5 and is seen to underestimate the lift appreciably at angles 
of attack greater than 3°. A somewhat better model of the flow at angle 
of attack is given by impact theory (ref. 4) which predicts increasing 
lift-curve slope with increasing angle of attack. This predicted varia- 
tion is plotted in figure 5; there is still a large gap between experiment 
and theory. 


The nonlinearities indicated in impact theory may be evaluated dif- 
ferently using the crossflow concept proposed in reference 5 and further 
developed in reference 6. If the crossflow contribution is added to the 
continuation of the initial lift-curve-slope line of figure 5, the sum 
agrees quite well with the experiment. This crossflow-lift increment was 
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calculated on the basis of ТО percent of the supercritical drag coeffi- 
cient of & circular cylinder being developed on an expanding body at 
angles of attack above 2.5. This value, ТО percent, was chosen on the 
basis of results given in reference 6 for the nose of an ogive cylinder. 
The agreement, perhaps better than might reasonably be expected, still 
indicates the possibility of accurate lift calculation by presently avail- 
able methods. | 


The effect of the spike on the lift appears to be small. Its frontal 
area, less than 4 percent of the model base area, would indicate that even 
if the lift carry-over from the nose were equal to the nose lift, the 
total lift would be increased by only à percent. The spike apparently did 
not affect the lift on the main body to any great degree, as was indicated 
by the success in calculating the total lift. 


The center-of-pressure position was estimated theoretically, ignoring 
the lift developed by the spike cylinder and assuming that the pressures 
on the aftercone were not affected by the presence of the cylinder. The 
resulting center of pressure, essentially the center of plan-form area of 
the conical portions, is predicted to be 74.2 percent of 2 aft of the 
nose. The median experimental value of 13.3 percent agrees well with this 
prediction, If the lift carry-over from the nose cone onto the cylinder 
is assumed equal to the nose-cone lift and acts at one third the cylinder 
length aft of the nose cone, the predicted location is 73.4 percent 1 
aft of the nose. It is expected that at large angles of attack, the center 
of pressure moves forward slightly since the spike-cylinder lift variation 
is expected to be more nonlinear than that for the conical portions. 


The estimates of the zero-lift drag of the projectile were based on 
references A, 7, and 8 for wave drag and reference 9 for skin friction. 
Ihe base pressure was &ssumed to be 10 percent of free-stream static 
pressure, based on extrapolation of the data of reference 10. The wave 
drag estimated by use of reference 4 was about 20 percent below that pre- 
dicted by use of reference 7 for а 5 cone. The value from reference 7 
only is presented in figure T. 


The theoretical drag build-up for both configurations is presented 
in figure T. The skin friction was calculated using local conditions and 
the effect of heat transfer was included. The data of reference 9, and 
subsequent data, ав yet unpublished, obtained by the ваше technique, were 
used in these calculations. 


The effect of heat transfer was to increase the skin friction by 
some 35 percent over that which would have existed without heat flow, 
according to reference 9. In order to see if the wall to free-stream 
temperature ratio was the same as for the tests of reference 9, estimates 
were made of the maximum temperature rise the model could have experienced. 
Before firing, the model temperature was around 5409 R and the free-stream 
static temperature was 300° R. At M = 8.6 the recovery temperature was 


жа 
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around 43009 R, and the flight duration was 0.008 second. At the outset, 
it is useful to note that even if the entire model surface were to be 
raised to the melting point, 14009 R, the temperature difference causing 
heat transfer would be reduced by only 23 percent. That this extreme 
temperature rise did not, in fact, occur is evidenced by the perfect con- 
dition of models observed at the last shadowgraph station. 


Calculations of the maximum possible temperature rise of the model 
based on the appendix of reference ll indicated a temperature rise on the 
spike cylinder, just behind the nose cone, of 309 R and 1000? R at the 
apex. Thus the temperature rise on 97 percent of the wetted area was 
less than 30° R and the maximum temperature rise at the nose, as indicated 
by its survival, was less than 880° R. Thus it is seen that the wall to 
stream temperature relation of the skin-friction tests also applies to 
the present test. 


Two calculations of the local Mach number and dynamic pressure, on 
which the skin-friction calculations were based, were made. The first 
was based on reference { and ignored the effect of the spike on the flow 
over the main body. The second was based on reference Т for the nose cone, 
and shock-expansion theory, reference 8, for the remainder of the body. 
The pressure distributions and skin-friction results from these two cal- 
culations were nearly identical, indicating no important effect of the 
cylinder on the pressure distribution on the remainder of the body at 
small incidence. 


Ав is indicated in figure T, boundary-layer transition generally 
occurred well downstream of the spike base. In many cases the flow was 
laminar over the entire body length, but it never remained so throughout 
an entire flight. Thus it appears that at the present test conditions 
of high Mach number and large heat-transfer rate into the model surface, 
& laminar boundary layer is extremely stable. 


CONCLUDING REMARKS 


Presently available theories may be used for engineering estimates 
of lift, center-of-pressure location, and wave drag of the projectile and 
the drag of the 5° half-apex-angle cone from which the projectile shape 
was derived. The lift for the projectile at angles of attack around 3° 
and above is noticeably greater than predicted by linear theory. 


A good estimate of friction drag was made on the basis of presently 
available date which include the effect of heat transfer. 


Ames Aeronautical Laboratory 
National Advisory Committee for Aeronautics 
Moffett Field, Calif., Aug. 23, 1954 
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Figure 3.- Photomicrographs of model surface and profile, X400. 
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(a) Angle-of-attack and vertical-position histories , test 3. 
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(b) Angle-of-attack history, test |. 


Figure 4.- Plots of basic data for two tests. 
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(c) 


Test 3, station 4 


(d) 


Test 4, station 1 


(e) 


Test 5, station 2 


Figure 6.- Shadowgraphs 
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(k) Test 11, station | 
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Figure 7.- Variation of drag of projectile and a 5? cone 


with position of transition location. Мж 8.6. 
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